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ABSTRACT 

Nearby clusters such as Virgo and Coma possess galaxy distributions which tend to be aligned with the principal 
axis of the cluster itself. This has also been confirmed by a recent statistical analysis of some 300 Abell clusters 
where the effect has been linked to the dynamical state of the cluster Moreover, the orbits of satellite galaxies in 
galactic systems like our own Milky Way also demonstrate a high degree of anisotropy - the so-called Holmberg 
effect, the origin of which has been the subject of debate for more than 30 years. This study presents the analysis 
of cosmological simulations focusing on the orbits of satellite galaxies within dark matter halos. The apocentres 
of the orbits of these satellites are preferentially found within a cone of opening angle ^40° around the major axis 
of the host halo, in accordance with the observed anisotropy found in galaxy clusters. We do, however, note that a 
link to the dynamical age of the cluster is not well established as both our oldest dark matter halos do show a clear 
anisotropy signal. Further analysis connects this distribution to the infall pattern of satellites along the filaments: 
the orbits are determined rather by the environment of the host halo than some "dynamical selection" during their 
life within the host's virial radius. 

Subject headings: galaxies: formation — cosmology: theory — methods: numerical 



L INTRODUCTION 

Observations of the distribution of bright elliptical galaxies 
within the Virgo galaxy cluster show a remarkable collinear ar- 
rangement (West & Blakeslee 2000). Moreover, this axis also 
appears to be part of a filamentary bridge connecting Virgo and 
the rich cluster Abell 1367. This phenomenon has akeady been 
recognized (Arp 1968; Bingegeli, Tammann & Sandage 1987) 
but only with accurate measurements of distance did it become 
possible to discriminate between this being a genuine three- 
dimensional structure or merely a chance alignment of galaxies. 
West & Blakeslee (2000) based their distances upon the surface 
brightness fluctuations method (cf. Blakeslee, Ajhar & Tonry 
1999 for a recent review), concluding that not only were the 
brightest cluster members distributed anisotropically over the 
cluster, but the distribution of dwarf ellipticals was elongated 
in the same direction (Bingelli 1999). Plionis et al. (2003) in- 
vestigated this substructure-cluster correlation statistically us- 
ing 303 Abell clusters. They affirm that there is indeed such a 
signal and also that this signal is related to the dynamical state 
(and the environment) of the cluster 

Moreover, even on smaller scales - in galactic systems - 
there are clear observational indications that the distribution 
of the orbits of satellite galaxies is biased towards the galac- 
tic pole (Holmberg 1969) giving an anisotropic distribution. 
Zaritsky et al. (1997) found that satellites of (isolated) disk 
galaxies are scattered asymmetrically about the parent galaxy 
and aligned preferentially with the disk minor axis. At first 
this signal was only prominent for distances out to ^50 kpc 
from the host (Holmberg 1969), but an analysis based on a 
much larger sample (cf. Zaritsky et al. 1993) of satellites ex- 
tends it to 200 kpc (Zaritsky et al. 1997). This result is con- 
firmed by a study of the satellites orbiting M3 1 (Hartwick 2000; 



Grebel, Kolatt & Brandner 1999). In addition, for the one 
galaxy where individual satellite orbits are known, the Milky 
Way, there is also evidence that that the orbits are preferen- 
tially polar (Zaritsky & Gonzalez 1999). This is derived from 
information based upon the alignment of satellites on the sky 
(Kunkel & Demers 1976; Lynden-Bell 1982), the orientation of 
the Magellanic Stream (Mathewson, Clearly & Murray 1974), 
the three-dimensional distribution of satellites (Majewski 1994; 
Hartwick 1996), and their actual velocities (Scholz & Irwin 
1994). 

All these observations clearly indicate that on scales span- 
ning from galaxy clusters down to galaxies there is a signal indi- 
cating a correlation between the alignment of substructure and 
the shape of the gravitational potential this substructure moves 
in. However, the source of this alignment of satellites has been 
debated for a number of years, with two potential solutions sug- 
gested to explain this puzzling arrangement. If the distribution 
of infalling satellites is initially spherical, then dynamical selec- 
tion may preferentially destroy or suppress the star-formation 
in those on more equatorial orbits^. For instance, Penarubbia, 
Kroupa & Boily (2002) have argued that in flattened dark mat- 
ter halos, the timescale for dynamical friction is substantially 
longer for satellites that are on polar orbits. Alternatively, the 
present non-isotropic distribution of satellite systems may re- 
flect the fact that it has always been non-isotropic, with the 
satellites being accreted along preferred directions. These hy- 
potheses can be tested in numerical simulations of structure 
formation. However, in previous attempts to decipher this sig- 
nal from cosmological A^-body simulations there still remains a 
certain amount of uncertainty (Zaritsky et al. 1997). The first 
fully self-consistent simulations targeting the subject were per- 
formed by Tormen (1997). They though lacked resolution in 



time, space and mass to perform a detailed analysis of the satel- 
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lite population for differing environments, and hence more re- 
fined simulations are required. Tormen (1997) could not fol- 
low the satelUte distribution within the host's virial radius but 
rather tracked all progenitors prior to accretion. This allowed 
him to analyze the infall pattern rather than the orbital evolu- 
tion of the satellites. The aim of this study is to investigate this 
subject with the latest state-of-the-art high-resolution A^-body- 
simulations. We focus upon a detailed analysis of the temporal 
and spatial properties of satellite galaxies residing within host 
galaxy clusters that formed fully self-consistently within a cos- 
mological framework. 

2. THE SIMULATIONS 

There is mounting evidence that, despite its problems, the 
Cold Dark Matter structure formation scenario provides the 
most accurate description of our Universe. Observations point 
towards the standard ACDM according to which the universe is 
comprised of about 28% dark matter, 68% dark energy, and lu- 
minous baryonic matter (i.e. galaxies, stars, gas, and dust) at a 
mere 4% level (cf. Spergel et al. 2003). This so-called "concor- 
dance model" induces hierarchical structure formation in which 
small objects form first and subsequendy merge to form larger 
objects. 

The analysis presented in this study is based upon a suite 
of eight high-resolution A'-body simulations within the ACDM 
concordance cosmology. They were carried out using the pub- 
licly available adaptive mesh refinement code MLAPM (Knebe, 
Green & Binney 2001). MLAPM reaches high force resolution 
by refining all high-density regions with an automated refine- 
ment algorithm. The refinements are recursive: the refined re- 
gions can also be refined, each subsequent refinement having 
cells that are half the size of the cells in the previous level. 
This creates an hierarchy of refinement meshes of different res- 
olutions covering regions of interest. The refinement is done 
cell-by-cell (individual cells can be refined or de-refined) and 
meshes are not constrained to have a rectangular (or any other) 
shape. The criterion for (de-)refining a cell is simply the num- 
ber of particles within that cell and a detailed study of the appro- 
priate choice for this number as well as more details about the 
particulars of the code can be found elsewhere (Knebe, Green & 
Binney, 2001). 

Each run focuses on the formation and evolution of one par- 
ticular dark matter halo containing of order greater one million 
particles. The force resolution reached by MLAPM in the high- 
density regions is Ih"^ kpc corresponding to roughly 0.25% of 
the host's virial radii /?vir- Following Lacey & Cole (1994) we 
define the formation redshift ziom as the time when the halo 
contains half of its mass today. The host halos were carefully 
chosen to sample a variety of different triaxiaUties and merger 
histories as summarized in Table 1. The virial radii given in 
that Table correspond to the point where the density of the host 
(measured in terms of the cosmological background density pb) 
drops below Ayir = 340 with Myir being the mass enclosed by 
that sphere. All satellite galaxies orbiting in and around the host 
halo are identified using the MLAPM-based halo finder described 
in a companion paper (Gill, Knebe & Gibson 2004). This new 
halo finder is based upon the MLAPM grid structure and works at 
the same resolution level as the A^-body code itself thus assur- 
ing not to miss or include any objects not resolved in the actual 
simulation. The satellites are then individually traced from Zform 
onwards until redshift z = 0. The mass spectrum of those satel- 
lites accounted for in the following analysis can be described 
by a declining power-law dn /dM oc M"" with a « 1 .7 - 1 .9 in 



the range from 2 x lO'"/;"' M© (applied mass-cut correspond- 
ing to 100 particles and explaining the rather 'low' number for 
A^sat(< '"vir) in Table 1) up to ~ lO^^/i"' Mq. The total fraction 
of mass locked up in satellites never exceeds <15% though. 

The high temporal sampling of the outputs (i.e. Af « 0.2 
Gyrs) allows us to accurately measure the orbital parameters 
of the satellite galaxies and a detailed analysis of their dis- 
ruption and survival history will be presented elsewhere (Gill, 
Knebe & Gibson 2004). Gill et al. show that the distribution 
of orbital eccentricities for the satellite population peaks about 
the value e w 0.6 where e = I- p/a is defined using the last 
pericenter p and apocenter a of its orbit. The pericenter distri- 
bution itself peaks about 35% of /?vir for all host halos and most 
of the satellites (i.e. more than 75%) had had one full orbit 
with the maximum number of orbits found to be 3-5 depending 
on the host halo. The angular momentum vector of the satel- 
lite ^^at'^^at' X ^sM° was calculated at its last apocenter passage 

'^sat ■ 

In order to probe the ahgnment with the shape and orientation 
of the host halo we calculated the eigenvectors £'1,2,3 (with Ei 
being the major axis) of its inertia tensor using only the "core" 
region as defined by the 6* refinement level in MLAPM, i.e. the 
boundary of this refinement level is an isodensity contour. Ac- 
cording to the refinement criterion adopted in the simulations 
the 6* level surrounds material about 3000 times denser than 
Ph or in other words 9 times denser than the material at the 
virial radius. The host halos are well described by the density 
profiles advocated by Navarro, Frenk & White (1997) with con- 
centration in the range c = 5 — 7. Therefore, a density of roughly 
9 X /5(rvir) corresponds to about the half-mass radius of the host. 
The eigenvectors now define the orientation of the host halo and 
the coordinate system used to measure the orbits of the satel- 
lites, respectively. Moreover, its eigenvalues a > b > c can be 
used to construct the triaxiality parameter T = (a^-b^)/(a^-c^) 
(Franx, lllingworth & Zeeuw 1991) presented in Table 1, too. 

The left panel of Fig. 1 presents the (normalized) distribution 
of angles between and £i,host- This graph shows that there 
is a clear trend in at least six of the eight halos for the two vec- 
tors to be aligned (the distribution peaks at 0° and 180°, respec- 
tively) meaning that the orbits of the satellites are preferentially 
found along the major axis of the host. The right panel of Fig. 1 
shows the (normalized) distribution of angles between the an- 
gular momentum vector of the satellite at its last apocenter 
and £1,^81. clearly revealing evidence for these two vectors to 
be parallel (distribution peaks at 90°). We Uke to stress that only 
satelHte galaxies that at least had one or more complete orbits 
were taken into account in Fig. 1; the figure is not based upon 
the infall pattern of satellites as investigated by, for instance, 
Tormen (1997). We also need to stress that the distributions 
presented in that Figure are normalized, i.e. they are corrected 
for the bias introduced by plotting them as a function of the 
angle 6 rather than cos0. The data is binned equally spaced in 
angle ranging from 0° to 180°, and therefore the area probed on 
the sphere varies with 6. Therefore the distribution needed to be 
normalized by the the respective area A specified by the actual 
range of angles [0- A0/2, 0+ A0/2]. This area is proportional 
to A oc cos(6l - A6»/2) - cos(6» + A6»/2). 

3. CONCLUSIONS AND DISCUSSION 

There are two possible scenarios that explain the 
substructure-cluster aUgnment found in observations of galaxy 
clusters (Plionis et al. 2003; West & Blakeslee 2000) and the 
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Table 1 

Summary of the eight host dark matter halos. 



Halo 


/?v.r[/!-' kpc] 


v^f^l km s-i] 


MvirLlO^V Mq] 


Zform 


age [Gyr] 


A^sat(<rvir) 


T 


# 1 


1349 


1125 


2.87 


1.16 


8.30 


158 


0.67 


#2 


1069 


894 


1.42 


0.96 


7.55 


63 


0.87 


#3 


1081 


875 


1.48 


0.87 


7.16 


87 


0.83 


#4 


980 


805 


1.10 


0.85 


7.07 


57 


0.77 


#5 


1356 


1119 


2.91 


0.65 


6.01 


175 


0.65 


#6 


1055 


833 


1.37 


0.65 


6.01 


85 


0.92 


#7 


1014 


800 


1.21 


0.43 


4.52 


59 


0.89 


#8 


1384 


1041 


3.08 


0.30 


3.42 


251 


0.90 



analysis of A^-body-simulations presented in this study. Either 
some dynamical process can be held responsible, meaning that 
initially the orbits were randomly distributed and only some 
satellites survive to give rise to the observed correlation in 
Fig. 1, or alternatively, their orbital parameters are imprinted 
upon them at the time they enter the host halo as already pointed 
out by Tormen (1997). 

"Dynamical destruction", whereby satellites on more equato- 
rial orbits suffer more tidal disruption, can be rejected because 
a comparable (although not as pronounced) alignment signal 
is found when restricting the analysis to those satellites that 
are in fact disrupted: We placed mass-less tracer particles at 
the last centre of the satellites before they were classified "dis- 
rupted". These "disrupted systems" are presented as the thin 
lines in Fig. 1. 

We therefore conclude that it must be the initial distribution 
of the satellite systems that is responsible for the present day 
alignment. For example, galaxies are "funneling along the fil- 
aments" (Kitzbichler & Saurer 2003) which has also been con- 
firmed by X-ray observations of the spatial distribution of sub- 
structure in galaxy clusters (West, Jones & Forman 1995) and 
biases the types of possible orbits. It would now be reassuring 
to confirm a link between the host halo shape and the surround- 
ing environment, i.e. are the filaments that feed the halo with 
material preferentially angled with respect to the orientation of 
the host halo? When projecting the positions of the satellites 
onto a sphere at the time they enter the virial radius of the host it 
is clear that they are not randomly distributed; they cluster in di- 
rections linked to the filamentary structure surrounding the host 
halo, as found before observationally (e.g. West & Blakeslee 
2003; PUonis & Basilikos 2002, and references therein) and in 
A^-body-simulations (e.g. Faltenbacher 2002; Hatton & Ninin 
2001; Onuora & Thomas 2000; Colberg et al. 1999; SpHn- 
ter et al. 1997; Tormen 1997). Tormen (1997) already pointed 
out that there is a strong alignment between the distribution 
of infalling satellite galaxies and the shape of the dark matter 
host. However, Tormen's simulations did not have the spatial 
and mass resolution to investigate the orbits of satellites within 
the host halo. We are, for the first time, tracing extremely well- 
resolved dynamics of the satellites within the dark matter hosts 
confirming that the alignment for both, disrupted and surviving, 
satellites is maintained for several orbital periods. Fig. 2 now 
presents a visual impression of the trajectories of all surviving 
satellites identified at formation time of host halo #3 until red- 
shift z = 0, qualitatively supporting the scenario that these orbits 



are linked to the filamentary structure which will leave its im- 
print in observed anisotropy. 

Plionis et al. (2003) also estabUshed a Unk between the align- 
ment and the dynamical age of the clusters. Our simulations 
though show that this finding can not be generalized: both our 
oldest halos, in which the satellites had as many as 4-5 orbits, 
do show the correlation signal. We do not observe a trend for a 
"randomization" of the orbits in older halos. The satellites ac- 
tually preserve the alignment with the host they had when they 
first fell into the cluster. 

In Section 1 we showed that the (observational) substructure- 
alignment signal spans from galaxy clusters down to galactic 
scales. Our results are mostly applicable to the former, and 
any extrapolation to smaller scales has to be handled with care. 
In the case for galaxies, for instance, the interaction between 
the galactic disk and the incoming satellite could enhance de- 
struction in the galaxy plane, i.e. satellite on prograde orbits 
decay faster than the ones on retrograde (or polar) orbits due to 
orbital resonances between the disk and the satellites (Peiiarub- 
bia, Kroupa & Boily 2002). Based on theoretical and numeri- 
cal studies (cf. Lacey & Cole 1993 and Moore et al. 1999) it 
should be possible to re-scale our simulations to a Milky Way 
sized object by requiring that the maximum of the circular ve- 
locity curves of our halos equals 220 km s~' . The scaling factor 
lies in the range of about 4-5 (cf. v^^^" values in Table 1) and 
hence our "rescaled dark matter halos" would correspond to 
the ones in the observational data with virial radii in the range 
260-295/!"' kpc. This scaling factor in length entails a scal- 
ing in mass of 64-125 (simply the length scale to the power 
of three). This brings the mass of our galaxy clusters down to 
~ 10 '^/j~' M0 which agrees with the dark matter mass inferred 
for our Milky Way (Freeman 1996). More problematic, how- 
ever, are the ages of our systems: our halos are only <8.3 Gyrs 
old opposed to 12 Gyrs for the Milky Way. Satelhtes in the 
Milky Way had the chance to complete nearly twice as many 
orbits leaving more space for an explanation based upon the 
dynamical destruction scenario. 

When re-scaling our data, putting aside the age issue, and try- 
ing to explain the Holmberg effect, another uncertainty comes 
into play: the orientation of the stellar disk. In a triaxial po- 
tential there are in general stable closed orbits about both the 
major and the minor axis (e.g., Binney & Tremaine 1987). So, 
in principal the disk plane could either be perpendicular to the 
major or perpendicular to the minor axis. Under the assumption 
it lies perpendicular to the major axis of the dark matter halo 
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the satellites in our simulations will be on polar orbits. There 
are indications that this configuration results in the most stable 
disk configuration within a triaxial halo (Hayashi et al. 2003). 
Moreover, even though there are clear indications that the an- 
gular momentum of the dark matter is well aligned with the 
minor axis of the halo (e.g. Warren et al. 1992), van den 
Bosch et al. (2002) showed that the angular momentum of 
the baryonic component (i.e. gas) not necessarily follows that 
of the dark matter distribution. They found an average mis- 
alignment between Lg;,s and Lqm of the order of 40° in their 
numerical simulations. Hence the orientation of the galactic 
disk with respects to the dark matter halo is not well deter- 
mined. Turning to observations does not resolve this ques- 
tion, with studies of polar rings indicating strongly oblate halos 



(i.e. lodice et al. 2003), while others suggest dark matter halos 
are more spherical, or even oblate (Oiling & Merrifield 2000; 
Ibata et al. 2001). The situation is not at all clear, but the as- 
sumption that disks are perpendicular to the major axis (which 
is in agreement with results presented by Hayashi et al. 2003) 
would provide an explanation for the Holmberg effect, even 
though this is a very speculative interpretation. 

The simulations presented in this paper were carried out on 
the Beowulf cluster at the Centre for Astrophysics & Supercom- 
puting, Swinburne University. The support of the Australian 
Research Council and the Swinburne Research Development 
Grants Scheme is gratefully acknowledged. GFL thanks Suede 
for "Introducing the band". 
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Fig. 1 . — The left panel shows the (normalized) distribution of angles between position vector of the satellites measured at the last apocentre and the major 
axis iJi.host of the halo, d values of 0° and 180° indicate alignment with the the principal axis of the host. The right panel shows the (normalized) distribution of 
angles between angular momentum vector L^Jj'" of satellites measured at the last apocentre and the major axis £i host of the halo. A 6 value of 90° means that the 
orbital plane of the satellite contains the major axis of the host where 9 of 0° and 180° would require for the plane to be perpendicular to the host's principal axis. 
The thin histograms are for disrupted sateUites. In both cases only sateUites more massive than 2 x 10"'/i"' Mq (100 particles) that had one full orbit were taken 
into account. 
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Fig. 2. — Orbits of all satellite galaxies identified at formation time (dark) to redshift 7. = (light). The left panel shows how the host halo fits into the surrounding 
large-scale structure presented at formation time. The right panel zooms into the region marked in the left panel, this time not showing low-resolution particles. 
The shape (and position) of the underlying host halo at redshift z = is indicated by the best-fit ellipse to the grid used to calculate its triaxiality. The vertical line 
extending to top edge of panel indicates the principal axis of the host. With the assumptions presented in this paper, the pole of the host galaxy within the dark 
matter halo is aligned with this principle axis. 
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